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Abstract 

The goal is to identify the class of distributions to which the distri- 
bution of the maximum of a Levy process with no negative jumps and 
negative mean (equivalently, the stationary distribution of the reflected 
process) belongs. An explicit new distributional identity is obtained 
for the case where the Levy process is an independent sum of a Brown- 
ian motion and a general subordinator (nondecreasing Levy process) in 
terms of a geometrically distributed sum of independent random vari- 
ables. This generalizes both the distributional form of the standard 
Pollaczeck-Khinchine formula for stationary workload distribution in 
the M/G/l queue and the exponential stationary distribution of a re- 
flected Brownian motion. 
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1 Introduction and preliminaries 

Let X = {X t \ t > 0} be a Levy process with no negative jumps. It is 
standard knowledge that in this case Ee~ aXt = e^^ 1 is finite for all 
a > and that 

tp(a) =ba + ^—a 2 + / (e~ ax - 1 + axl {x < 1} ) v{dx) (1) 

2 «/(0,oo 
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where b is real, a is nonnegative and v is a measure satisfying 

(x 2 A l)u(dx) < oo (2) 



L 



(0,00) 

where a Ab = min(a, b). 

It is also well known that —EX t /t — <fi'(0) = — b + , xv(dx) 
and in particular EX t is well defined and can be either finite or +00 
but never —00. In particular if EX t < (equivalently <p'(0) > 0) then 
f(i 00) xv {dx) < b so that in particular b > and ^ xv(dx) < 00. 
Therefore, in this case ip has the following form 

2 p 

<p(a)= fia + ^-a 2 + (e~ ax - 1 + ax) v{dx) (3) 

where fi — <p'(0) = b — J, x ^ xv{dx) > 0. 

When in addition ^ xv{dx) < 00 then the Levy process is an 
independent sum of a Brownian motion and a pure jump subordinator 
(nondccrcasing Levy process). In this case, ip becomes 

2 f 

ip(a) = ca + %-a 2 - (l - e^) u(dx) , (4) 

2 J (0,00) 

where c = fj, + J. Q ^ xv{dx). 

It is also well known that for any Levy process with no negative 
jumps and ip' (0) > 0, if we denote M = sup t>0 X t , then 

Ee-" = ^ (5) 

and that the limiting and stationary distributions of the (Markov) 
process W t — X t + L t , where L t = — info< s <t X s , is the distribution of 
M. We (and others) call the formula the generalized Pollaczek- 

Khinchine (PK) formula. The reason for the name is that for the 
special case of the M/G/l queue, the underlying driving (Levy) process 
is a compound Poisson process (a subordinator) minus t for which the 
Laplace-Stieltjes transform of the limiting and stationary distribution 
of the workload process is the celebrated Pollaczek-Khinchine formula 
and may be found in virtually all basic queuing theory textbooks. 

There are various textbooks where Levy processes are discussed 
and where the above results may be either found directly or concluded 
from (e.g., [1, 4, 8, 13] and pages 19-34 of [12]). There are quite a few 
different proofs in the literature for the generalized PK formula, mostly 
via the application of the Wiener-Hopf factorization (e.g., [6, 15]), weak 
convergence (e.g., [5, 14]) and martingales (e.g., [7]), but this is not 
the scope here. 
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For some recent work on the distribution (rather than Laplace- 
Stieltjes transform) of a-stable Levy processes see [2, 9, 10] and further 
references therein. For the case with phase type upward jumps (and 
general negative jumps) see [11]. For some other results see also [3]. 
These papers also include an extensive list of references to texts and 
further related literature. 

One is not required to be a Levy process expert to read this paper 
and all the knowledge which is needed for what follows is covered above. 
In particular, this may easily be taught in any course where Levy 
processes are touched upon. 

2 The case of an independent sum of a 
Brownian motion and a subordinator 

Let X be a Levy process which is an independent sum of a Brownian 
motion and a subordinator having a finite mean. In this case, as seen 
in the previous section, the (Laplace-Stieltjes) exponent can be written 
in the form 

2 2 p 

¥>(<*) =ca + - / (1 - e-°*Mdx) , (6) 

1 J(0,oo) 

where D = J, Q . xv(dx) = J °° v{x, oo)dx < c. Denoting 

F e (a)= r e -™«^ds, (7) 
Jo v 

p = ^ < 1 and A = , we observe that the exponent may be rewritten 
like this 

<p(a) = ca (l + j - P F e (aj) (8) 

and in particular (p'(0) = c(l — p) so that the generalized PK formula 
has the form 

<V(Q) = i-p f9 x 

<p{a) l + f- P F e (a) [) 

If a 2 = (A = oo) and v(0, oo) < oo then one obtains the well known 
(original) PK formula for the M/G/l queue. When v(0, oo) = oo it is 
interesting to observe that exactly the same formula is valid without 
change only that F e is then the stationary excess life distribution as- 
sociated with the jumps of the subordinator and unlike in the renewal 
process setting, its density approaches oo in the neighborhood of zero. 
If the Levy measure is zero then p — and the formula becomes the 
Laplace-Stieltjes transform (LST) of an exponential distribution with 
rate A. This is also well known to be the LST of the stationary distri- 
bution associated with a one dimensional reflected Brownian motion 
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with negative drift. The interesting discovery is that these two results 
can be unified. To see this, note that 



A l-p 

(10) 

n— v 7 

and so we have the following result. 

Theorem 1 For a Levy process with no negative jumps satisfying 
f( 1 j xv(dx) < oo, (p'(0) > and with the notations defined above, 
let N ~ G(l - p) in the sense that P[N = n] = (1 - p)p n . Let 
Xq , Xi , X 2 , . . . ~ cxp( A) f= for A = oo ) and Y\ , Y 2 , . . . ~ F e and as- 
sume that N, X a ,Xi, . . ., Yi,Y 2 , . . . are all independent. Then °J^^ 
is the LST of the following random variable 

N 

X + J2(Xn+Y n ) (11) 

n=l 

where an empty sum is zero. 

When v i then we are left with X as expected (the Brownian motion 
case) and when a 2 I we are left with X^=i a l so as expected 
(distributional form of the PK formula). 

An interesting special case occurs when the jumps of the Levy pro- 
cess have a phase-type distribution. Then the residual life also has a 
phase-type distribution and thus also X-i + Yi and Xo + ^^li^n + ^n)- 
This a special case of the results reported in [11] but with an easier 
proof but also a far less complicated setup. 

It is very easy to check that also a converse holds. 

Theorem 2 Assume that < p < 1, 0<A< oo, / is nonnegative, 
nonincreasing with J °° f(y)dy = 1 (possibly with f(x) — > oo as x 10) 
and denote F(x) = J* f(y)dy. Let N <~ G(p), Xq,X\,... ~ exp(A) 
(= for X = oo) and Y\,Y 2 , ■ ■ ■ ~ F where all random variables are 
independent. Then there exists a Levy process which is an independent 
sum of a Brownian motion and a subordinator having a negative mean 
for which ^1 ls the LST of 

N 

X + J2(Xn + Y n ) , (12) 

n=l 

where an empty sum is zero. This Levy process is unique up to a 
constant time scale. 



<V(0) 
ip(a) 
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Proof: To see this we can simply perform reverse engineering. First 
assume without loss of generality that / is right continuous, otherwise 
we take its right continuous version which gives the same F. Let 
v((x, oo)) = Bf(x) for any constant 8 > and thus v((a, b]) — 8(f(a) — 
f(b), which uniquely characterizes v. Note that 



xv(dx) = I v((x,oo))dx = I 6f(x)dx — j3 < oo (13) 

(0,oo) JO JO 

as required. Recalling that p — ^, that necessarily p = 1 — p and 
since v = 8 we must set c = Now, if A = oo we set a 2 = and 

otherwise, from A = ^§ we set a 2 = X (i- P ) ■ Finally, 

ip'(0) = c-u=j^- p>Q (14) 

so that all requirements are met. 

The fact that this Levy process is unique up to a constant time 
scale is evident since if 



<V(o) _ ojy/(o) 

Lp(a) V(a) 



(15) 



then necessarily ip{a) ~ 795(a) for 7 = , so that if X is a Levy 
process with exponent <p and Y with tp, then {X lt \ t > 0} is distributed 



like Y. 



3 The general spectrally positive case 

The fact that ^ggl i s the LST of some proper distribution for any 
Levy process with no negative jumps and a negative mean is deduced 
indirectly from its various proofs which always exploit the connection 
with either the supremum or the reflected process. However, we are 
not aware of a more explicit and direct way of showing this than The- 
orems 1,2 together with the following. 

Observe that if J^ Q ^ xv(dx) = 00 then the Levy process is not a 
sum of a Brownian motion and a subordinator. In this case (since 
(p'(0) > and thus f, t ^ xv{dx) < 00), we recall from (3) in Section 1 
that 

tp(a) = pa + + / (e- ax - 1 + ax) v{dx) (16) 

2 Ao.oo) 

and that (p'(0) = p. Now, with 

22 r 

<p e {a)=p J a+^-+ (e~ ax - 1 + ax) u{dx) (17) 

2 J(e,oo) 
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for e > 0, wc clearly have that ip e (a) — > (p(a) as e 4- 0, that (p' e (0) = 
if' (0) = yu and thus 

Um ^ = VW (18) 

and that aip , ^ — > 1 as a i 0. 

Now (/? e is the exponent of a Levy process which is an independent 
sum of a Brownian motion and a a compound Poisson process for which 
we have identified the distribution associated with in Theorem 1. 

From this the following is immediate. 

Theorem 3 A distribution has the generalized PK EST for some Levy 
process with no negative jumps and a negative mean if and only if it 
belongs to the closure of the family of distributions defined in Theo- 
rem 2. 



References 

[1] Applebaum, D. (2009). Levy Processes and Stochastic Calculus, 

2nd ed., Cambridge University Press. 
[2] Bcrnyk, V., Dalang, R. C. and Pcskir G. (2008). The law of the 

supremum of a stable Levy process with no negative jumps. Ann. 

Probab. 36, 1777-1789. 

[3] Chaumont, L. (2010). On the law of the supremum of Levy pro- 
cesses. arXiv:1011.4151vl [math.PR]. 

[4] Bertoin, J. (1996). Levy Processes, Cambridge University Press. 

[5] Czystolowski, M. and W. Szczotka. (2010). Queueing approxima- 
tion of suprema of spectrally positive Levy process. Queueing Sys- 
tems 64, 305-323. 

[6] Harrison, J. M. (1977). The supremum distribution of a Levy 
process with no negative jumps. Adv. Appl. Probab. 9, 417-422. 

[7] Kella, O. and W. Whitt. (1992). Useful martingales for stochastic 
storage processes with Levy input. J. Appl. Probab. 29, 396-403. 

[8] Kyprianou, A. E. (2006). Introductory Lectures on Fluctuations 
of Levy Processes with Applications, Springer. 

[9] Kwasnicki, M., Malecki, J. and M. Ryznar. (2011). Suprema of 
Levy processes. arXiv:1103.0935vl [math.PR]. 

[10] Michna, Z. (2011). Formula for the supremum distribution of a 
spectrally positive a-stable Levy process. Stat. Probab. Letters 
81, 231-235. 

[11] Mordecki, E. (2002). The distribution of the maximum of the Levy 
process with positive jumps of phase-type. Theory Stock. Process. 
8, 309-316. 



6 



[12] Protter, P. E. (2003). Stochastic Integration and Differential 
Equations, 2nd ed., Springer. 

[13] Sato, K. I. (1999). Levy Processes and Infinitely Divisible Distri- 
butions, Cambridge University Press. 

[14] Szczotka, W. and W. A. Woyczyriski. (2003). Distributions of 
suprema of Levy processes via heavy traffic invariance principle. 
Probab. Math. Stat. 23, 251-272. 

[15] Zolotarev V.M. (1964) The first passage time to a level and the be- 
havior at infinity of processes with independent increments. The- 
ory Probab. Appl. 9, 653- 661. 



7 



